QCD Corrections to QED Vacuum Polarization 



P. Castelo Ferreira 
pedro . castelo . f erreiraOist .utl.pt 

J. Dias de Deus 
jddOf isica. ist .utl .pt 

CENTRA, Instituto Superior Tecnico, Av. Rovisco Pais, 1049-001 Lisboa, Portugal 



Abstract 

We compute QCD corrections to QED calculations for vacuum polarization 
in background magnetic fields. Formally, the diagram for virtual ee loops is 
identical to the one for virtual qq loops. However due to confinement, or to the 
growth of as as decreases, a direct calculation of the diagram is not allowed. 
At large we consider the virtual qq diagram, in the intermediate region we 
discuss the role of the contribution of quark condensates (qq) and at the low- 
energy limit we consider the tt", as well as charged pion tt+tt^ loops. Although 
these effects seem to be out of the measurement accuracy of photon-photon 
laboratory experiments they may be relevant for 7-ray burst propagation. In 
particular, for emissions from the center of the galaxy (8.5 kpc), we show that 
the mixing between the neutral pseudo-scalar pion ttq and photons renders a 
deviation from the power-law spectrum in the TeV range. As for scalar quark 
condensates (qq) and virtual qq loops are relevant only for very high radiation 
density ~ 300 MeV/ fw? and very strong magnetic fields of order ~ 10^"* T. 
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1 Introduction 



In the presence of background electromagnetic fields second order QED correc- 
tions in the finc-stnicture constant a to vaciiuni polarization due to quantum 
vacuum oscillations, i.e. electron-positron virtual loops [1,2], is a well stud- 
ied subject, including Delbriick Scattering [3], photon-splitting [4,5], photon- 
photon interactions [6, 7] and to semi-classical interactions with pseudo-scalar 
particles [8-10]. QED effects are well established, Dclbruck scattering has been 
experimentally observed in light scattering by heavy nuclei [11] as well as con- 
tributing in second order in aZ to the Lamb shift [12]. However, a direct 
signature for interaction with axion-like pseudo-scalars has been consistently 
verified [13,14]. Independent experiments to detect this sort of interaction have 
now been considered [15,16]. In addition in astro-particle observations, due 
to the wide range of energies accessible, QED effects are the main contribu- 
tion for the optical depth. Such effects include particle production (also known 
as photon desintegration) [17], vacuum polarization and photon splitting [18]. 
Also astro-particle observations seem to be the most promising candidate to test 
pseudo-scalar interactions [19-21] trough, for example, analysis of gamma-ray 
burst conversion rates [22] and its polarization characteristics [23] . 

In this work we study the second order perturbative corrections in the fine- 
structure constant due to the strong interactions. Namely we analyse quark 
loops, quark condensates and meson contributions to vacuum polarization. The 
second order contribution to the polarization of the vacuum due to electron- 
positron virtual pair production is given by the Euler-Heisenberg Lagrangian [1, 
2] 

4 {F^.F'^-^f + 7 (e'^-^'PF^^FspY^ 
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The respective vacuum polarization dispersion relation for radiation propagating 

in vacuum under an orthogonal magnetic field Bq -C B'^ is expressed in terms of 
the eigenvalues and A jj , respectively orthogonal and parallel to the magnetic 
field [4,5], 

a;^,ll = k(l-Xl„Bi) , 

(2) 

Al = S^.Bl , A| = Ui.Bl . 

Also within the framework of QED we have the contribution of other fermionic 
loops than the ones due to the electron. Relevant to the present study we have 
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the muon loops fj,fi which give the contribution 



AC^ = f=[:^] = 5.43 X IO- 



CS) 



being of the same order of magnitude of the QCD corrections that we are ad- 
dressing here. Due to its higher mass the r gives a contribution five orders of 
magnitude lower ('--^ 6 x 10^^^), hence not being relevant here. 



2 QCD Contributions 

Naively, we can expect that the same kind of physics applies to quark-antiquark 
virtual pair production. In the presence of an external field we have in general 
the diagram of figure [T] We write in the case of qq virtual pair production in 





Figure 1: The diagrams for fermion-antifermion loops and the exchange of a vr' 
neutral meson. The vertex n^^y includes the axial anomaly. 



order the Euler-Heisenberg Lagrangian 
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The factor of 3 comes from the summation over colours and Qq stands for the 
quark fractional charge. As for the quark masses vriq correspond to the renormal- 
ized masses that appear in the quark propagator. Under an external magnetic 
field Bq the relation of the polarization due to qq with the polarization due to 
ee corresponding to electron-positron loops is given in terms of the corrections 
to the parallel and orthogonal vacuum dispersion relation eigenvalues and 
A| by 



(5) 



Here we consider the up quark mass « 5MeV and charge Qq = 2/3 and 
Sq = WAq/wtot < 1 is a phase space correction due to confinement of strong 
interactions. We know that at low energies there are no free quarks, therefore 
quark loops carrying small momenta cannot be considered. The way out is to 
introduce a lower cut-off A, in the loop momenta such that only the higher 
momenta contribution to the loop is considered. 

The probability for the full range of momenta (i.e. g]0, +oo[) is given by 
the series [24] 

wtot = —Y- Y] . (6 

n=l 

Due to confinement and the increase of as for small values of , we introduce 
a cut-off Aq that truncates the series ^ by excluding the low region 

aB^ \ - 1 _ g 



obtaining the relation 



-A,= (^) • (8) 



For the light quarks (summing over up and down quark masses) with mass of 
order m, ~ lOMeV [25] we have riA, 3600 such that Sq ~ 10-i°'*/l^l. This 
value is obtained from the leading term contribution from the above series for 
o^A, • Here we considered A, ~ 600 MeV, this being the value for which the 
strong interactions coupling constant becomes unity, as 1 [26], such that 
below this energy threshold, QCD is in a non-perturbative regime. The free 
quark loop contribution to vacuum polarization is therefore negligible. This 
contribution will only be relevant for very strong magnetic fields of order B ^ 

The low-energy quark states (corresponding to the light mesons) are the tt's. 
In low energy physics these particles can be used as fundamental bosons within 
the framework of Chiral Perturbation Theory (ChPT) [27]. Therefore below the 
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cut-off < hq the main contribution is due to the neutral meson tt*^ with an 
effective scalar Lagrangian [2] : 

(9) 

= = 2.49 X 10-2 GeV-^ , 

where the coupling gT^-yy is taken from the Adler-Bell-Jackiw anomaly coeffi- 
cient [28] and the pion decay constant is taken to be /^r = 93AIeV. The re- 
spective contribution to the dispersion relation of radiation traveling in vacuum 
(corresponding to the tt diagram of figure [1]) is 

Xf - , Xf - A^,oA| , 

2 A (10) 

We considered the pion mass m.^ ~ 135 MeV, for higher masses the contribu- 
tions are for most applications negligible when compared to the tt'^ effect. 

In addition we can have loops of the lighter charged mesons n~^TT~ whose 
Euler-Heisenberg Lagrangian reads [2] 

(11) 

45(B,-)2 ' ^ eh ' 

hence contributing a correction to the vacuum dispersion relation of approxi- 
mately the same order of magnitude given by 

Af = ^A^.i Al , Af* = ^A^.i Af 

(12) 

There is yet another contribution that we can consider. In the presence of 
background magnetic fields there is a vacuum polarization contribution due to 
quark condensates. Within the Schwinger- Euler-Heisenberg formalism [1,2] in 
the context of ChPT [29] a correction to the vacuum dispersion is obtained 

Al = ACcAl , A| = , 

= 8e: = i^^n^j (13) 
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Taking the quark condensate ultra-violet cut-ofF A « 300 MeV we obtain that 
A^c ~ 1-69 X 10~^°. Next we give some details on how quark condensates are 
obtained and explain which regimes exist depending on the loop momentum. 
The parallel vacuum polarization for ChPT is given by the integral 



n 



f 
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ds I, 
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a B cot(a B s) 
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(14) 
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This distribution is represented in figure [51 The contributions considered here 




Figure 2: (a) The integrand The poles at s = {n — l)n/aB (for n = 1,2, . . . ,oo) 

are marked by vertical lines and contribute to the pion vacuum polarization, (b) The 
same integrand between both cut-offs = 135 MeV and A = 300 MeV for B = 5.5 T . 
It corresponds to the marked region between the poles at s = and s = n/aB of (a). 



are due to the poles below the cut-ofF s < 1/A^. For weak fields the only 
pole that contributes for pion loops is at s = 0. It corresponds to the tt+tt" 
loops and the relative magnitude of its effect has already been discussed and 
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Table 1: The several QCD effects in the presence of weak fields and their mag- 
nitude for the several ranges of the loop momenta p. To exist, the quark con- 
densate requires a high density of energy. 

is given in equation p^ . Above the cut-off s > we consider the quark 

loops instead of the meson distributions. The novel interesting feature in this 
framework is that we have a new contribution between the pole s — Q and 
s = 7r/a|i?|, it corresponds to the quark condensate. We note that from a 
more fundamental level based in Nambu-Jona-Lasinio theory [30] the quark 
condensates contribution is of the same order of magnitude [31]. There is an 
important point to stress, NJL consider explicit actions for the quarks instead of 
the effective actions for the mesons considered in ChPT, the condensate cut-off 
A should correspond in NJL to the confinement energy. These theories were 
originally motivated by superconductivity and the relation between ChPT and 
NJL is equivalent to the relation between Landau-Ginzburg effective theory [32] 
and the Bardeen-Cooper-Schriffer microscopic theory [33] for superconductivity. 

We summarize the allowed effects and their magnitude for several ranges of 
loop momenta p in table [TJ 

We note that the ttq is a ^ pseudo-scalar such that CP-symmetry is 
conserved, the Lagrangian ([5]) is a scalar. More generally we may consider also 
the contributions of other ^ pseudo-scalars as the r/'s, however their masses 
are higher than the mass of the pion (the lighter being to^ w 547 MeV) , hence 
their contribution is negligible by several orders of magnitude. 

As for quark condensates are O"*"*" scalars. Although we have already present 
the effects of these condensates in equation (jl3p obtained within the mean- 
field framework of ChPt, we note that we can recast this effects diagrammat- 
ically considering the scattering of photons by an intermediate scalar. How- 
ever if one demands CP- invariance the effective action is expressed as £ = 
9{qq)4'{qq)P"^"'^ Pp.i'y wherc (f>(qq) rcprescnts now a scalar. In the same manner 
one may consider the other 0++ scalar mesons /o and oq [34]. Although these 
processes are allowed and the couplings are of the same magnitude of the pion, 
are negligible due to much higher masses [25]. 

In addition we also note that the quark condensates may only exist when 
very high densities of energy are present (E) ^ 300MeV/ fni^ [29] . These val- 
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ues are only accessible in very dense plasmas (for example in neutron stars [35]) 
or for very high fluxes of high energy radiation. For example for radiation 
energies in the TeV range propagating in vacuum are required fluxes over 
10^^ photons ■m~'^ . It is understood that these fluxes include both the prop- 
agating and the background photons. 

3 Vacuum Birefringence 

The relevant radiative corrections that contribute to vacuum birefringence are 
of second order in the fine-structure constant [4,5]. The usual classical wave 
equation in order is linear in the photon field A [5]. Hence for a static 
transverse magnetic field Bq, gathering the results from the previous section we 
obtain, due to QCD effects, a correction on the refractive indexes eigenvalues 
given by [4,5,15] 

X± 

\\ 

The directions || and ± correspond respectively to the parallel and transverse 
directions to the external magnetic field, is given in ^ and the several 
corrections are ordered by magnitude significance according to the estimatives 
given in ([3]), pil)) . ([T^ and (|13p . The above equations result in having different 
refractive indexes for the parallel and perpendicular directions to the magnetic 
field [15] 

iV|| = l-iA|| , iV^ = l-iA^, (16) 

which introduce a phase shift in the propagating wave. Considering a linearly 
polarized wave of wave number \i — koz which polarization makes an angle 
00 with a static magnetic field Bq both gain an cUipticity and its polarization 
is rotated due to the vacuum effects (see for example [15]). The polarization 
rotation is given by 

A6I = i (A|| - A_l) Az sin(26io) , (17) 

being Az the distance traveled by the radiation, and the cUipticity is given by 

V' = -^(A|| -Ai)Azsin(20o) , (18) 

being lu the radiation frequency. The relative magnitudes of rotation induced in 
vacuum by the several effects presented here are pictured in figured We note 
that the contribution to the rotation of the pseudo-scalars have the same sign 
of the fermionic loops contributions (A|| > Aj^), while the contribution due to 

the charged pion loops have the opposite sign (Aj|* < X^)- 
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1 + AC^ + A^e + Ia^^± ] 6 Bo' , 



= 14 l + Ae^ + Ae.o + -A^,±UeB, 



(15) 



Figure 3: Relative contributions to the polarization rotation as given by equation p7| ) 
from the several effects as a function o/A^i, i.e. each effect magnitude in relation to the 
magnitude effect due to electron-positron loops (ee). Both axes are in logarithmic scale. 
The continuous line coincides approximately to the PVLA experimental conditions [13] 
with Bo = 5.5 T and Az = 10^ m. The marked points are labeled and correspond to 
the QED corrections due to electron-positron loops (ee), the muon-antimuon loops 
(up.) interchange of the neutral pion (tv^), quark condensates ({qq)) and charged pion 
loop (n'^n~ ). The ellipticities are obtained by re-scaling these results by the radiation 
frequency to. 



We note that the usual QCD scale is set by Kqcd ~ 2QQMeV, however 
as already discussed, we also know that for energies of approximately Ag = 
600 MeV the strong running coupling constant as is of order of unity and the 
perturbative regime of QCD is no longer valid [26] . Therefore the correct value of 
the cut-off corresponding to low-energy quark condensate is not exactly known 
and should be in the range 200 < A < 600 MeV . This value should correspond 
to the chiral phase transition energy of the Nambu- Jona-Lasinio theory [30] . 

As already put forward in [9] (see also [10]) these effects, in particular the 
contribution of the pion, is several orders of magnitude lower than the one due to 
the electron vacuum oscillation as well as of the theoretical axion contribution. 
In addition we recall that the quark condensate contribution is only present if 
very high density of radiation is considered. 

We conclude that vacuum birefringence due to QCD corrections is negligi- 
ble for all known physical systems. As already put forward before in [9], the 
main contribution is due to the neutral pion being many orders of magnitude 
below the vacuum polarization effects of virtual electron loops which, although 
a well established phenomena within QED, has not directly been detected any 
polarization rotation, neither in laboratory experiments, neither in astrophys- 
ical environments. However in the next section we give an example where the 
light scattering by the 7r° meson may have measurable effects. 
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4 High Energy 7-ray Propagation 



We can also apply the results derived so far to the propagation of 7-ray burst 
in background magnetic fields where the the effects studied in the previous 
sections seem to be relevant due to the very high radiation energy. The main 
contribution to the attenuation of the 7-ray spectrum for high energies (E above 
10^^ eV) is photon desintegration (particle-antiparticle pair production) due 
to interaction with the background electromagnetic fields [17]. These effects 
result in a exponential decaying law for the 7-ray spectrum, ~ E~^'^^\ For 
radiation from the center of the galaxy (corresponding to z = 8.5 kpc) the 
value of the decaying exponent is F « 2.25 [22]. The radiation flux is of order 
10~^ photons TO~^ TeV~^ [19] (corresponding to a radiation density of order 
p.y ~ IQ"^^ MeV/ fm^). Therefore for this particular case, the only relevant 
contribution discussed in this work is from the 7r° meson which we address 
next. 

We note that it is also expected that the axion-like pseudo-scalar contri- 
bution has visible effects in the high energy range (of order of TeV) either 
increasing or decreasing the optical dept depending on the values of the mass 
and photon coupling constant considered. More specifically it is expected that 
for more stable particles with long decaying time (low decaying rate) as the light 
axion, the optical depth increases [9,20,21], while for less stable particles with 
lower decaying time (greater decaying rate), the optical depth decreases. What 
distinguishes between these cases is the relation of the pseudo-scalar mass (m^) 
to the photon pseudo-scalar coupling [g-y^,) properly taking in consideration the 
background and traveling radiation energy. We will return to this discussion by 
the end of this section. Specifically the equation for a generic pseudo-scalar cf) 
mixing with the photon is [9, 20, 21] 



{uj - id:, + M) 



A 



A, 



= 



(19) 



with 



M 



^11 + All 
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and the several entries given by 
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(20) 
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The approximation in A-^^ corresponds to a linearization of the cross section 
for photon desintegration in the TeV range of the energies (E) for distances 
Zo « 8.5 kpc and we have considered the photon mass negligible (m^ « 0). 
The approximation in S\\ ± corresponds to neglecting the contributions from /x/x 
loops and the charged pion Tr~^TT~ loops. As derived in [21], for non-polarized 
radiation in gaussian magnetic field distributions, the conversion probability of 
photons to pseudo-scalars is 

1 / 3Po^ \ 

^7-0 = 3 (1 - e ) , 

f X 2 (22) 

Po « 0.4xlO-M^^^j . 

These expressions are applicable to the mixing of the neutral pseudo-scalars to 
photons in 7-ray bursts with coupling given by g = g^^/lO"^ Ge^"^, the root 
mean square magnetic field strength Bq = 1 ^Gauss, the energy Ew given in 
units of \QTeV and the mass given in MeV. z is the distance to the source in 
pc and s the size of the magnetic field domains also in pc [21]. For 7-ray burst 
from the center of the galaxy one has z — 8.5 kpc and s = O.Olpc [21, 22]. For 
the particular cases of the ttq mixing we have g — 2.49 x 10'* and m ~ 135 Af eV . 
The deviation to the power law considering this effect is pictured in figure ID 




Figure 4: Deviations from the power law due to ttq exchange. The open circles and 
filled circles are data points from HESS collaboration corresponding to the data sets 
from July/august of 2003 and 2004 [^^J (^^^ o,lso [21]). The filled line represents the 
power law best fit dN/dE ~ E~^ with T = 2.25 and the dashed line to the contribution 
of the TTQ . 



As we have put forward in the beginning of this section, depending on 
the specific coupling constants, masses and existing energies the effects of the 
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pseudo-scalar mixing to photons renders quite different results. The relevant 
expressions to compare correspond to the diagonal and off-diagonal compo- 
nents of the matrix for pseudo-scalar photon mixing. Respectively the relevant 
expressions to compare are rn^^zj s and g^BE. This analysis was originally car- 
ried in [9] and also considered in [20,21]. The main differences between this 
several works are the numerical values for the parameters of the propagation 
equation. While in [9] the length z/s traveled by radiation is of order of Km, 
in [21] are considered astrophysical environments traveling several kpc which 
allows to saturate the scalar-photon oscillations as given by equation . The 
quantitative similarities of our results in relation to [21] is simply due to the 
ratios of the couplings to mass squared of the heavy axion considered there and 
the pion considered here, being different only by about one order of magnitude 
21 X g^/mj « 9 axion /"m axion- fo'" [20] we note that it is considered a very 
light axion m'^^^g^ ^ g-y axionEB such that the massless limit niaxion ^ is 
taken. This last case is clearly not applicable to the pion where we have that 
> g^T^BE and explains why for a relatively heavy intermediate pseudo- 
scalar (with low decaying times) the optical-depth is decreased while for a very 
light intermediate pseudo-scalar (with high decaying times) the optical-depth 
in increased. 

It could also be considered the more generic case of several pseudo-scalar 
and scalar mixing. Not only the theoretically suggested axion, as well as the 
quark condensate effects when are met the conditions for their formation (for 
example close to neutron stars and magnetars [35]). 

5 Conclusions 

In this work we have computed the QCD corrections to QED vacuum polariza- 
tion effects. Although we conclude that the contribution to vacuum birefrin- 
gence of the effects presented here are negligible when compared to the effect of 
virtual electron loops they have observable consequences for high energy 7-ray 
propagation. In particular we have shown that the neutral pion mixing with 
photons significantly contribute to a deviation from the power-law spectrum 
in the TeV range that may be relevant when considering the superposition of 
other pseudo-scalar effects in this range [19,21]. As for quark condensates and 
virtual quark loops we deduced that only for very high radiation energy fluxes 
(Pc > 300Mey//TO3) and strong magnetic fields {B > 10"T), their effects 
may be relevant. Hence near neutron stars and magnetars [35] these effects may 
affect 7-ray polarization [23] . 
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